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Abstract. Using a full BABAR data sample of 426 fb we present improved measurements of the ratio 
7?(D'*') = S{B -> D^'h-Vr)! S(B -> D'-'H-Vf), where I is either electron or muon. We measure HiD) = 
0.440 ± 0.058 ± 0.042 and HiD') = 0.332 ± 0.024 + 0.018. These ratios exceed the Standard Model predictions 
by 2. Oct and 2.1a, respectively. The results disagree with the Standard Model predictions at the level of 3.4it. 
The ratios are sensitive to new physics contributions in the form of a charged Higgs boson. However, the access 
cannot be explained by a charged Higgs boson in the type II two-Higgs-doublet model. 



1 INTRODUCTION 

The semileptonic physics in B meson sector played a 
prominent role in investigating of new physics effect at 
low-energy region. Semileptonic transitions are the sim- 
plest process in B mesons decay. In the Standard Model 
(SM), the heavy b quark decays to either a c or an m 
quark and the virtual W boson |l-3]. Experimentally, 
semileptonic decays have the advantage of large branch- 
ing fraction and are used to determine the weak cou- 
plings, the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements, \Vcb\ and \Vub\ iRl- 

The parton level diagram of B — > decays 
where D**^ refers to either a D or a D* meson is shown 
in Fig.[T] The decay B £)**V is sensitive to charged 
Higgs contribution at the tree level. The three-body de- 
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Figure 1. The parton level diagram for B — > D**'t Vr decays 

cay B — > permits the study of decay distribution 

which discriminate between and H exchange jsl 0] . 

The decay of B — > D^*'t Vr with r lepton in the final 
state offer possibilities of significant new physics contribu- 



tions that is not present in the process where light lepton 
such as electron and muon in the final state. The study of 
B — > have already shown the new physics con- 

tributions can be over-constrained ifrl-fiol. The existing 
studies of the B —> D**V Vt based on two-Higgs-doublet 
model (2HDM) predict a substantial impact on the ratio 
and «(Z)) 10 HE [11. 
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2 The BABAR DETECTOR AND DATA SET 

The BABAR detector was operated at the PEP-II 
asymmetric-energy storage rings at the SLAC National 
Accelerator Laboratory. The data used in this analysis 
were collected with the BABAR detector. We analyze 
data recorded with the BABAR detector at a center of 
mass energy of 10.58 GeV. The data sample consist of an 
integrated luminosity of 426 fb^\ corresponding to 471 
xlO^BB pairs. An additional sample of 40 fb taken at 
energy 40 MeV below the T(4S) resonance. This addi- 
tional sample of data is used to study the continuum back- 
ground from the decays of e^e' — » qqiy) pairs where q 
can be u, d, s, c, r. 

A detail description of the BABAR detector is pre- 
sented elsewhere Ol ill . The momenta of the charged par- 
ticles are measured in a tracking system consisting of a 
5-layer double sided silicon vertex tracker (SVT) and a 
40-layer drift chamber (DCH). The SVT and DCH oper- 
ate within a 1 .5 T solenoid field and have a combined solid 
angle coverage in the center of mass frame of 90.5%. A 
detector of internally reflected Cerenkov radiation (DIRC) 
is used for charged particle identifications of pions, kaons, 
and protons with likelihood ratios calculated from dE/dx 
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measurements in the SVT and DCH. Photons and long- 
hved neutral hadrons are detected and their energies are 
measured in a CsI(Tl) electromagnetic calorimeter (EMC). 
For electrons, energy lost due to bremsstrahlung is recov- 
ered from deposits in the EMC. 

3 ANALYSIS 

In this analysis, instead of measuring the absolute branch- 
ing fraction of B — > D^*h^Vr, we measure the ratios 



<B{B D^*H-ve) 



(1) 



where ( is either electron or muon. In the standard model 
(SM), the relative rate HiD^*^) have less than 6% uncer- 
tainty llT2i1 . In the decay of B — > D***t v^, we construct the 
T lepton only from the purely lepton decays: e^fieVr 
and e^v^Vj so that the signal events {B D^*W^Vj) 
and the normalization events (B — » D^*^{^V() are identified 
by the same particles in the final state. When taking the ra- 
tio of 'R{D^*^), the various sources of uncertainties will be 
canceled and reduced. 

We reconstruct candidate events produced in T(45 ) — > 
BB decays by selecting the hadron decay of one of the B 
meson {B,ag). The other candidate events are reconstructed 
semileptonically, specially a charm meson (either charged 
or neutral D or D*) and a charged lepton (either e or yu). 
The signal events {B — > Vt) and the normaliza- 

tion events {B — > D^*H~ve) are extracted using unbinned 
maximum-likelihood fit to the two-dimensional distribu- 
tions of the invariant mass of the undetected particles. Ba- 
sically it is the invariant mass of the neutrinos. 



= Pmiss = iPe-e- ' PB,„, ' Pd'-> ' Pf) 



(2) 



where p is the four-momenta of the colliding beams, B,ag, 
D^*' and charged lepton, respectively. The lepton three- 
momentum in the B rest frame is denoted by p*^. The 
distribution of the lepton three-momentum of the signal 
events is softer than the distribution of the lepton three- 
momentum of the normalization events because the ob- 
served lepton in the signal events is a secondary particle 
originated from the t decay, — > {^V(Vt- 

If all particles are properly reconstructed the invariant 
mass of the undetected particles {m^ . ) with a single miss- 
ing neutrino peaks at zero, whereas the signal events which 
have three missing neutrinos have a wide inl^.^^ distribution 
that extends from -1 GeV^ to 10 GeV^. The two observ- 
able kinematic variables are used to select the Btag candi- 
dates: 

' (3) 
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= M - 

Y beam 



Ptag 



and 



AE = E, 



lag 



(4) 



where the p,ag and E,ag refer to the center-of-mass momen- 
tum and energy of the B,ag. Eheam is the center-of-mass of 
a single beam particle. If the B decays are correctly re- 
constructed, the distribution of the niEs is centered at the 



B meson mass with a resolution of 2.5 MeV. The distribu- 
tion of AE is centered at zero with a resolution of 18 MeV. 
In this analysis we required niEs > 5.27 GeV and \AE\ < 
0.072 GeV. 

The main background contributions to the signal 
events are the following: 

• The decay of B — > D**(t^ /£')v( where D** mesons refer 
to the charm resonances heavier than the D* meson such 
as Dq, Di, D'j, and orbital excitations of the cq pairs. 
The decay of B — > D**('V( where the D* meson decays 
to D^*^7fi peaks in the m^^^.^^ distribution. These events 
are estimated using the Monte Carlo samples. 

• Charge cross-feed events: these events come from the 
decay of B — > D**(t )vf. These background events 
were reconstructed with the incorrect charge where one 
of the charges particles in the final state has been as- 
signed to the different B meson. 

• Other BB background: These events come from the de- 
cay of B ^ £)(t,**)£)K«*)+ j-jjg large leptonic and 
semileptonic branching fractions of Z)^ mesons. We es- 
timate these events using Monte Carlo sample and its 
contribution is fixed in the fitting process. 

• Continuum background: these events come from the de- 
cay of e^e^ —> qqiy) pairs where q can be u, d, s, c, t. 
To estimate this background we use the additional data 
sample of 40 taken at energy 40 MeV below the 
T(45' ) resonance. 

Figure |2] shows the distributions of w^^.^^ and p*^ variables 
for the signal events (6%), normalization events (82%), 
B — » D**{v (4%), and combinatoric and continuum back- 
ground (8%). All distributions are normalized to 1000 en- 
tries after all selections are applied. 
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Figure 2. The comparison of the data and the fit projections 
for the four samples: the signal events (6%), normalization 
events (82%), B D"{v (4%), and combinatoric and contin- 
uum background (8%). The iues is on the left and the is on 
the right. All distributions are normalized to 1000 entries after 
all selections are applied. 



The results of the signal and normaUzation yields are 
extracted using unbinned maximum-likelihood fit to the 
two dimensional, mf^.^^-p*^ contributions. The fit is per- 
formed simultaneously to the four D^*^{ samples and four 
D^'^tt"^ samples. The distributions of each D^*^{ and 
D^'^T^i sample is fitted to the sum of either eight or six 
contributions, respectively. The fit relies on8x4-t-6x4 = 
56 probability density functions (PDFs). The two dimen- 
sional, mi l -p*f contributions for each of the 56 PDFs are 
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described in detail using smooth non-parametric kernel es- 
timator 013 1. 

The m^^j.^j distributions of the signal events and the nor- 
malization events cane be easily distinguished due to the 
different number of neutrino as the undetected particles in 
its corresponding decays in the final state. However, the 
"^miss distributions of the backgrounds resemble those of 
the signal events, and therefore in the fitting procedure 
these contributions are either fixed fitted or constrained 
by the D***7r"f Monte Carlo sample. Figure |3] shows the 
yield of B — > D*tv and the comparison of the /«,^„ „ and p*^ 
distributions of the B D*tv (data points) with the pro- 
jections of the results of the isospin-unconstrained fit. The 
region above the dashed line is the background component 
corresponds to BB background. 

The region below the dashed line corresponds to the 
continuum background. Due to the charged cross-feed 
events the . distribution peak around . = in the 

miss ^ miss 

background component. However, in the p*^ distributions, 
we only include events with mL- > 1 GeV^. 



Table 1. The yield results for the B D'tv channel where the 
uncertainties are the statistical only. 



Mode 




D'tv 


no*) 

S(D*Tv) 


639 + 62 245 + 27 
0.32 + 0.03 0.36 + 0.04 
1.71+0.17 1.74 + 0.19 


888 ±63 
0.33 ± 0.02 
1.76 ±0.13 


Table 2. The yield results for the B — » Dtv channel where the 
uncertainties are the statistical only. 


Mode 


D"tv D+tv 


Dtv 


^ sig 

S(Dtv) 


314 ±60 177 ±31 
0.43 ± 0.08 0.47 ± 0.08 
0.99 ±0.19 1.01 ±0.18 


489 ± 63 
0.44 ± 0.06 
1.02 ±0.13 



below the dashed line corresponds to the continuum back- 
ground. In the p*^ distributions, we only include events 
with ■ > 1 GeV^. Table |2] shows the fit results of the 
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Figure 3. The comparison of the m^„i^, (left) and (right) dis- 
tributions of the B — » D'tv (data points) with the projections 
of the results of the isospin-unconstrained fit. The region above 
the dashed line is the background component corresponds to BB 
background. The region below the dashed line corresponds to the 
continuum background. In the distributions we only include 
events with «J,^„j, > 1 GeV-. 



Figure 4. The comparison of the in^^j^^ (left) and (right) dis- 
tributions of the B — > Dtv (data points), with the projections 
of the results of the isospin-unconstrained fit. The region above 
the dashed line is the background component corresponds to BB 
background. The region below the dashed line corresponds to the 
continuum background. In the p'^ distributions we only include 
events with > 1 GeV^. 



We extract the branching fraction ratios as define in the 
following 

^(Z)(*)) ^ ^ ^ _ffi^ (5) 

^ norm ^norm 

where Nsig and Nnom are the number of signal and nor- 
malization events extracting from the fitting process, re- 
spectively. The esig/e„orm is the ratio of the efficiencies 
of the signal and the normalization events. We impose 
the isospin relations of ??(D*) = 'R(D*^) = ^(0*°) and 
HiD) = 9?(D+) = ??(D"). Table [U shows the fit results of 
the yield of B — > D*tv with the statistical uncertainties 
only. Figure |4] shows the yield of Z? — > Dtv compari- 
son of the OT^„^5 (left) and p*^ (right) distributions of the 
B — > Dtv (data points) with the projections of the results 
of the isospin-unconstrained fit. 

The region above the dashed line is the background 
component corresponds to BB background. The region 



yield of B ^ Dtv with the statistical uncertainties only. 

4 SYSTEMATICS UNCERTAINTIES 

The largest systematic uncertainties in this analysis is due 
to the poorly understood of the decay B — » D**({/t)v 
background. The systematic uncertainty due the PDF that 
describe these contributions including the uncertainty on 
the branching fractions of the four B —> D**(v decays, the 
branching fraction ratio of B ^ D**tv to B —> D**(v, and 
its relative efficiency. We assign 2.1% on fiiD) and 1.8% 
on HiD*), respectively. 

We also assign a systematic uncertainty due to the 
observed variation of the decay of B ^ D*ri(v, non- 
resonance B ^ D*n{n)£v, B D"(£/t)v, and D" 
D^*^nn. They are 2.1% for KiD) and 2.6% for KiD"). 

The other largest systematic uncertainties are due to 
the continuum and BB backgrounds. We assign 4.9% for 
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'R{D) and 2.7% for 'RiD*). The systematic uncertainties 
due to the PDFs for the signal and normalization decays 
ai-e 4.3% for •R(D) and 2.1% for The systematic 

uncertainties due to the efficiency ratios €sig and e„o„„ are 
2.6% on HiD) and 1.6% on n{D*), respectively. 

By choosing the decay of t lepton only from the purely 
lepton decays: — > e^v^Vj and — > e^v^v^, uncertain- 
ties due to the particle identification, final state radiation, 
soft-pion reconstruction, and other related detector perfor- 
mance are largely cancel in taking the ratios. They only 
contribute about 1 % in the systematic uncertainty. 

All systematic uncertainties are added in quadrature to 
assign the total systematic uncertainty. There is a positive 
correlation between some of the systematic uncertainties 
HiD) and 'RiD*). As a result the correlation of the total 
uncertainties is reduced to -0.27 for HiD) and 'RiD*). 

5 CONCLUSIONS 

We have measured the 'RiD) = 0.440 + 0.058 + 0.042 
and 'RiD*) = 0.332 + 0.024 + 0.018. These ratios ex- 
ceed the Standard Model predictions by 2. Oct and 2.7cr, 
respectively. They are disagree with the SM prediction 
at the level of 3.4cr. The results are compatible with the 
results measured by the Belle Collaboration [14, 15). To- 
gether with the results measured by the Belle Collabora- 
tion, it could be an indication of new physics processes in 
B mesons decay. 

The measured values of 'RiD) and 'RiD*) match the 
predictions of the particular Higgs model where the ra- 
tio of neutral Higgs field vacuum expectation values, tan/3, 
and the mass of the physical charged Higgs boson, 
is tan/3/t7iH+ - 0.44 ± 0.02 GeV"' and tan/S/mn* - 
0.75 ± 0.04 GeV"', respectively. Figure |5] shows the com- 
parison results of this paper with the charged Higgs boson 
of type II 2HDM predictions. The Standard Model predic- 
tions correspond to tan/S/niH* = 0. However, these results 




are not compatible with a charged Higgs boson in the type 
II 2HDM with a 99.8% confidence level for any value of 
tan/3/ mH+. More general charged-Higgs models or other 
New Physics contributions may explain the access of 'RiD) 
and 'RiD*). 
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Figure 5. The comparison results of this paper (blue band) with 
predictions that include a charged Higgs boson of type II 2HDM 
(red band). The SM corresponds to tan/iImM* = 0. 



